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ABSTRACT 

Context. The broad emission features in the Fe-Ko- region of X-ray binary spectra represent an invaluable probe to constrain the 
geometry and the physics of these systems. Several Low Mass X-ray binary systems (LMXBs) containing a neutron star (NS) show 
broad emission features between 6 and 7 keV and most of them are now interpreted as reflection features from the inner part of an 
accretion disk in analogy to those observed in the spectra of X-ray binary systems containing a Black Hole candidate. 
Aims. The NS LMXB GX 349+2 was observed by the XMM-Newton satellite which allows, thanks to its high effective area and good 
spectral resolution between 6 and 7 keV, a detailed spectroscopic study of the Fe-Ko- region. 

Methods. We study the XMM data in the 0.7-10 keV energy band. The continuum emission is modelled by a blackbody component 
plus a multicolored disk blackbody. A very intense emission line at 1 keV, three broad emission features at 2.63, 3.32, 3.9 keV 
and a broader emission feature in the Fe-Ka region are present in the residuals. The broad emission features above 2 keV can be 
equivalently well fitted with Gaussian profiles or relativistic smeared lines (diskline in XSPEC). The Fe-Ko- feature is better fitted 
using a diskline component at 6.76 keV or two diskline components at 6.7 and 6.97 keV, respectively 

Results. The emission features are interpreted as resonant transitions of S xvi, Ar xviii, Ca xix, and highly ionized iron. Modelling 
the line profiles with relativistic smeared lines, we find that the reflecting plasma is located at less than 40 km from the NS, a value 
compatible with the inner radius of the accretion disk inferred from the multicolored disk blackbody component (24 ± 7 km). The 
inclination angle of GX 349+2 is between 40° and 47°, the emissivity index of the primary emission is between -2.4 and -2, and the 
reflecting plasma extends up to (2-8) xlO** cm. 

Conclusions. We compare our results with the twin source Sco X-1 and with the other NS LMXBs showing broad relativistic lines 
in their spectra. We conclude that the blackbody component in the spectrum is the primary emission that hits the inner accretion disk 
producing the emission lines broadened by relativistic and Doppler effects dominant around the neutron star. 

Key words, line: identification - line: formation - stars: individual (GX 349+2) — X-rays: binaries — X-rays: general 



^^O 1. Introduction in atoll sources and from the horizontal branch (hereafter HB) 

, X . ^ , to the normal branch (NB) and to the flaring branch (FB) in Z 

^\ Low-mass X-ray bmaries (LMXBs) consist of a low-mass star 

(M < IMq) and a neutron star (NS), which generally has a rel- 



sources. 



GX 349+2, also known as Sco X-2, was called an odd-ball 



atively weak magnetic field (B < 10 G). In these systems, the , ^ ' r— — '- — — -i k^^^ i ^. ., 

V • ju ^- f a .-u among the Z sources flKuuIkers & van der Klia 19980 . Similar 

X-ray source is powered by accretion of mass overflowing the , ° r ^ , — ^ ^ — 1 , 

r> u 1 u e^u • » J f • J- 1 to the case of Sco X-1, GX 349+2 shows a short and under- 

Roche lobe of the companion star and forming an accretion disk , , r^, 

. ^, , , r A^vr. . • • . . /ivTc developed HB (if at all). The source variability in the frequency 

around the neutron star. LMXBs containing a neutron star (NS . , , , ■ . ■ ^ . . ■' 



d T \/tvD \ n A- ■ A A ■ ^ ^ A ^^ a range below 100 Hz is closely Correlated with the source position 

'■^ LMXBs) are generally divided into Z and Atoll sources, accord- ° . , „ ^ ■ . ,. .,, 

... »u »u J u • V 1 1 t-.- on the X-ray CD, as in other Z sources. Quasi periodic oscilla- 

ing to the path they describe in an X-ray Color-Color Diagram . . „ ,, ^„^x , ^ ,. , xtt^ x- 

/r^T^\ u A • . A- Itt — '■ D J T^i- I tions at kHz frequencies (kHz QPO) were detected in the NB of 

(CD) or hardness-intensity diagram OHasinger & van der Klisl .„ o, o,-. . i^^o. .. 

\ ]rioA ui J u • r~\ A i 1 Its Z-track Zhang, Strohmayer, & Swank 1998). However, GX 

119891) assembled by using the source count rate and colors cal- ^ ,. , - . , ' , , . 

, . ■ . • 1 / 11 T OA 1 V 349+2, which sometimes shows broad noise components chang- 

culated over a typical (usually 2-20 keV) X-ray energy range. . ' ... . . . , ^ , , r ■ r 

. ^ „ „ , » • J u 1 1 1 1 ins not only with the position in the Z, but also as a function of 

Atoll sources are usually characterized by relatively low lu- , ° . ■ , , , . ... ,.rr 

ni n o r \ j u . • /u u • the position in the hardness-intensity diagram, ditters somewhat 

minosities (0.01-0.2 Lfj,/ ) and often show transient behavior, „ , ,r, .. ,,, 

, -1 ry ,.u /-^ 1 *u . 1 • from the other Z sources and shows similarities to the behav- 

while the Z sources in the Galaxy are among the most lumi- . • , . , „ , ,o , . -, , 

T i\/rvT> ■ » ti 1 » T-AA- . 1- tor seen in bright atoll sources, such as GX 13+1 and GX 3+1 

nous LMXBs, persistently accreting close to the Eddmgton limit rrz — ^---^ — ■ — i N^^r. i i ^,xT ... nrrrrTm u^^^ k 

n :K Mc f ■ A- A ^ JKuulkers & van der Kid 1998tKyNeill et al. 200ll 2002h . 

(i^Edd) for a 1.4 Mq NS. The position of an individual source ^ — . . . — =r' ^' " *' ^ 

in the CD, which determines most of the observed spectral and | Di Salvo etal. | (|200lD, using BeppoSAX data, showed that 

temporal properties of the source, is thoug ht to be an indica tor of *e source energy spectrum below 30 keV could be wefl fit 

the instantan eous mass accretion rate (e.g.ll fasinger et alJ [T990t by a blackbody (with a temperature of 0.5-0.6 keV) and a 

Ivan der Khsl[T995l for a review). It has been suggested that the Comptonized component (with seed-photon temperature of 1 

mass accretion rate (but not necessarily the X-ray luminosity) keV and electron temperature of 2.7 keV). Three discrete fea- 

of individual sources increases along the track from the top left tures were observed in the spectrum: an emission line at 1 .2 ke V, 

to the bottom right, i.e. from the islands to the banana branch probably associated to L-shell Fe xxiv or Ly-a Ne x, an emission 

line at 6.7 keV (Fe xxv) and an absorption edge at 8.5 keV, both 

Send ofprint requests to: R. laiia, e-maiV. iariaOfisica.unipa. it corresponding to emission from the K-shell of highly-ionized 
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Fig. 1. Epic-pn lightcurve of GX 349+2 in the energy band 1-10 
keV. The bin time is 30 s. 



iron (Fe xxv). Ilaria et all ( |2004|) . analysing a long BeppoSAX 
observation, found similar parameters of the the continuum com- 
ponents below 30 keV, detected the presence of emission lines 
associated to Fe xxiv and Fe xxv, and detected an emission line 
at 2.65 keV associated to S xvi. The long BeppoSAX observa- 
tion allowed to study changes of the parameters of the Fe xxv 
emission line along the position of the source in the CD from 
the NB/FB apex to the FB inferring that the equivalent width of 
the line decrease from 77 eV to 18 eV going from the NB/FB 
apex to the FB. The width of the Fe xxv emission line, modelled 
by a Gaussia n line, was between 250 and 300 eV. 

Recently ICackett et alj (l2008b . analysing two Suzaku obser- 
vations of GX 349+2, modelled the Fe xxv emission line with 
a relativistic smeared line (diskline in XSPEC). The authors 
constrained the energy of the line to be between 6.4 and 6.97 
keV finding that the line energy was 6.97_o.o2 keV, associated 
to Fe xxvi and not compatible with a Fe xxv iron line a s previ - 
ously obtained bv Ilaria et all ( |2004 and lDi Salvo et al.1 (lioOlh . 
The emissivity index, yS, was -4.1+0.3, the inner radius in gravi- 
tational radii, Rg, was 8.0+0.4, the inclination angle of the source 
was 23 + 1 degrees, the equivalent width of the line was 76 + 6 
eV. 

In this paper we analyse an XMM-Newton observation of 
GX 349+2 finding that the prominent relativistic line is well con- 
strained around 6.7 keV and identified as the resonant transition 
of Fe xxv. In the following we discuss our results and compare 
them with the recent literature. 

2. Observation 



The XMM-Newton Observatory dJansen et all 1200 Ih includes 
three 1500 cm^ X-ray telescopes each with an European Photon 
Imaging Camera (EPIC, 0.1-15 keV) at the foc us. Two of 
the E PIC imaging spectrometers use MPS CCDs (JTurner et al.l 
1200 lb and one uses pn CCDs ( Striider et al.ll200ll). Reflection 
Grating Spectrometers (RGS, 0.35-2.5 keV. Iden Herder et all 
[2001 ) are located behind two of the telescopes. The region of sky 
containing GX 349+2 was observed by XMM-Newton between 
2008 March 19 16:42:41 UT to March 19 22:58:55 UT (OBSid 
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Fig. 2. Color-Intensity diagram of GX 349+2. The [4- 
10keV]/[l-4keV] hardness ratio increases for a count rate lower 
than 1580 c/s. The bin time is 100 s. 



05061 10101) for a duration of 22.5 ks. Since GX 349+2 is one 
of the brightest X-ray sources, in order to minimize the effects 
of telemetry saturation and pile-up the MOS 1 and MOS2 instru- 
ments were switched off" and the EPIC-pn camera was operated 
in Timing mode with medium filter during the observation. In 
this mode only one central CCD is read out with a time reso- 
lution of 0.03 ms. This provides a one dimensional (4'. 4 wide) 
image of the source with the second spatial dimension being re- 
placed by timing information. The faster CCD readout results in 
a much higher count rate capability of 1500 cts/s before charge 
pile-up become a serious problem for point-like sources. The 
EPIC-pn telemetry limit is approximatively 450 c/s for the tim- 
ing modeQ. If the rate is higher, then the counting mode is trig- 
gered and part of the science data are lost. This is the case of 
the observation presented in this work; saturation occurs for a 
55% of the observing time, giving an exposure time of 22.5 ks 
for RGSl and RGS2, respectively, and 10 ks for EPIC-pn. 

We plot the Epic-pn lightcurve in the 1-10 keV energy range 
in Fig. [1] which clearly shows that the count rate is largely vari- 
able. To understand whether the variability of the source inten- 
sity is energy dependent we extract the lightcurves in the 1-4 
keV and 4-10 keV energy range, respectively, and plot the hard- 
ness ratio [4-10 keV]/[l-4keV] versus the count rate in the 1-10 
keV energy range in Fig.|2] We find that the spectrum becomes 
harder for a count rate below 1580 c/s corresponding to the time 
interval between 10"* and 1.5 x 10"^ s from the start of the obser- 
vation. Plotting the 1-4 keV, 4-10 keV lightcurve and the [4-10 
keV]/[l-4keV] hardness ratio versus time (see Fig. |3]l, it seems 
that there is no correlation between the 1-4 keV and 4-10 keV 
count rate. The 1-4 keV count rate is quite constant at 1200 c/s 
between and 9 x lO"' s and between 1.5 x lO'* and 2.2 x lO"* s 
from the start, and it decreases between 9 x lO-' and 1 .5 x 10'* s, 
showing some flaring activity. The 4-10 keV count rate is quite 
constant at 530-550 c/s between and 6 x lO-' s and between 



' see Tab. 3 in XMM-Newton Users Handbook published on 15 July 
2008 
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Fig. 3. Upper Panel: 1-4 keV Epic-pn lightcurve. Middle Panel: 
4-10 keV Epic-pn lightcurve. Bottom panel: [4-10keV]/[l- 
4keV] hardness ratio. The bin time is 30 s for each lightcurve. 
The red dot-dashed vertical lines indicate the three selected time 
intervals. 



1.9 X 10"* and 2.25 x 10"* s, while increases to 600 c/s between 
6x10^ andl.9xl04 s. 

Since the count rate of the source is variable, we divide the 
lightcurve into three time intervals as indicated in Fig.[3]by the 
red dot-dashed lines. We base our choice on the 4-10 keV count 
rate that is 550 c/s in the first and third time interval and it is 
around 600 c/s in the second time interval. Our choice allows 
to highlight possible changes in the Fe-Ko' region of the spec- 
trum. In Tab. [1] we report the exposure times of the selected time 
intervals for RGSl, RGS2 and EPIC-pn instruments. 
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Fig. 4. Patterns of single (blue color) and double (red color) 
events extracted from the source region excluding RAWX 37 and 
38. The patterns are plotted versus the 60-2000 ADU channel 
range, corresponding to 0.3-10 keV. 



the source position (RAWX= 37). Background events were ob- 
tained from a box of the same width with RAWX between 2 and 
18. 

The 0.5-12 keV Epic-pn count rate extracted from the source 
region is almost 1 800 c/s, since the maximum Epic-pn count rate 
should be 800 c/s to avoid a deteriorated response due to photon 
pile-ufB we expect presence of pile-up. We checked the pile- 
up fraction at different energies using the epatplot tool in SAS 
finding a photon pile-up percentage of 1.3%, 2%, 3.2%, 5.1%, 



Table 2. Percentage of photon pile-up at different energies 



Table 1. Time intervals selected in the EPIC-pn lightcurve 





Time interval 


RGSl 


EPIC-pn 




s - TSTART[s] 


ks 


ks 


INT. 1 


0-6000 


6.0 


2.6 


INT. 2 


6000-19000 


13.0 


5.8 


INT. 3 


> 19000 


3.5 


1.6 



Note — The observation was divided into three time intervals as shown 
in Fig. [3] In the first column we report the temporal boundaries of each 
interval with respect to the start time of the observation. In the sec- 
ond and third column we report the exposure times, in kiloseconds, for 
RGSl, RGS2, and EPIC-pn during the three time intervals. 



2.1. Spectral analysis of the averaged spectrum 

We extract the X-ray data products of RGS and EPIC-pn cam- 
era using the Science Analysis Software (SAS) version 8.0.0 
extracting only single and double events (patterns to 4) from 
EPIC-pn data. Initially we extracted source EPIC-pn events from 
a 69.7" wide column (RAWX between 29 and 45) centered on 



Energy 
(keV) 


Pile-up (%) 
RAWX [29:45] 


Pile-up (%) 
RAWX [29:36] and RAWX[39:45] 


5 


1.3 ±0.1 


0.2 ±0.2 


6 


2.0 ±0.1 


0.6 ±0.2 


7 


3.2 ±0.2 


1.6 ±0.2 


8 


5.1 ±0.3 


1.9 ±0.4 


9 


7.4 ± 0.4 


4.8 ±0.5 


10 


9.8 ±0.6 


6.5 ±0.8 



Note — The associated errors are reported at 1 cr. In the first col- 
umn we report the energies, in the second, and third column the es- 
timated percentage of photon pile-up extracting the events from the 
whole source region and excluding the brightest columns at RAWX= 37 
and RAWX= 38, respectively. 



7.4%, and 9.8% at 5, 6, 7, 8, 9, and 10 keV, respectively. The 
values with the associated errors are reported in Tab. |2] Since 
our aim is the study of the Fe-Ka region of the spectrum we 
investigated how to minimize the photon pile-up effects in the 
energy range 5-8 keV. Initially we extracted the source EPIC- 
pn spectrum excluding the brightest column (RAWX= 37) in 

- see Tab. 3 in XMM-Newton Users Handbook published on 15 luly 
2008 



4 



R. laria et al.: A ionized reflecting skin in GX 349+2 



Table 3. Averaged Continuum and Identified Broad Gaussian 
Lines below 6 keV. 



Parameters 


Single Fe-Kci- 
Em. Line 


Fe xxv-Fe xxvi 
Blending 


N„ (10^2 cm-2) 


" ■^ -^-0.003 


736+° °°' 

"•'■'"-0.004 


klBB (keV) 
Nfifl (10-') 


1 052+''°i'' 
— u.uuy 


^•'°'-0.014 

l-057+°°lJ 

— U.U14 


^Tdiskbb (keV) 

^DISKBB 


1 05+°°^ 
210+"' 


1.04 ±0.03 
220 ± 20 


Eifev (keV) 
(Titrv (eV) 

hk,v (10-' phcm-2 s-') 
Equiv. width (eV) 
significance (cr) 


1 051+"°°* 

25 ±3 

22+:^ 
14 


1 05 1+°°°'* 
85+? 

-6 

24!^ 
21 ±2 
18 


E5 XVI (keV) 
Cs XVI (eV) 

Isxv/(10-3phcm-2 s-') 
Equiv. width (eV) 
significance (cr) 


2.62 (fixed) 
140 (fixed) 

9 -1+0.9 

Z,.J„ J 

6±2 
7 


2.62 (fixed) 

140 (fixed) 

2 5+°-^ 
"^■-'-0.6 

6±2 
7 


Ea^xv/;/ (keV) 

CAr XVIII (eV) 

lArXv;;/ (10-^ ph cm-- s-') 

Eqtiiv. width (eV) 

significance (cr) 


T 0^+0.03 

190+™ 

3 7+' -' 
-'•'-0.4 
11-1-^ 

IZ + J 

15 


T 09 +0.03 

200!^° 

4 0+' " 
^■"-1.1 

1/1+4 

6 


Ecn XIX (keV) 
(Tca XIX (eV) 

Ic«x/x(10-' ph cm-2 s-i) 
Equiv. width (eV) 
significance (cr) 


Q QQ+0.03 

ioo!l« 
9 9+0.7 

^■"^-04 

9±3 
9 


3.93 ± 0.02 
110+'° 

9 Q+(J?^ 

^■■J-0.4 
9+3 

9 


Xijd.o.f.) 


1.51(269) 


1.51(268) 



tween 2.1 and 2.5 keV from our analysis. The EPIC-pn spec- 
trum has been rebinned not to oversample the energy resolution 
by more than a factor 4 and to have 20 counts per energy chan- 
nel. The adopted Epic-pn energy range is 0.7-2.1 keV and 2.5-10 
keV. The spectrum is fitted using XSPEC version 12.3.1 

The extracted RGSl spectrum does not show instrumental 
issues and we select data in the 0.6-2 keV energy band for our 
analysis. The extracted RGS2 spectrum shows calibration prob- 
lems in CCD 9 and we will not use data from this instrument for 
our analysis. Since we are interested in the study of the Fe-Ka 
region, in the following we concentrate out spectral analysis on 
the Epic-pn data, after having checked that the RGS 1 spectrum 
was consistent with the Epic-pn data and does not show more 
features. 

The continuum emission is fitted adopting the same model 
recently proposed by Cackett et al. (2008) for the Suzaku spec- 
trum of GX 349+2: i.e. a blackbody component plus a multi- 
colored disk blackbody (DISKBB in XSPEC) both absorbed by 
neutral neutral matter (WABS in XSPEC); because of the nar- 
rower energy range of our Epic-pn data (0.7-10 keV) with re- 
spect to the Suzaku data we did not add to the model a power- 
law component necessary to fit the Suzaku data above 10 keV. 
We obtained a reduced x^g^/ - 7.73 with 282 degree of free- 
dom (d.o.f.). The residuals, plotted in Fig. |5] suggest the pres- 
ence of broad emission lines at 1, 2.7, 3.35, 3.9 keV and, finally 
a broader emission line between 6 and 7 keV. Initially, we added 




T3 
0) 



Note — Best-fit values of the parameters obtained for the aver- 
aged spectrum fitting the emission features with Gaussian profiles. 
Uncertainties are at the 90% confidence level for a single parameter. The 
parameters are defined as in XSPEC. The continuum emission is fitted 
with a blackbody plus a multicolored disk blackbody, diskbb. The 
four broad emission features are modelled adopting Gaussian compo- 
nents. All components are subject to interstellar absorption. 



o 



the CCD finding that the photon pile-up fraction did not change 
significantly; we obtained a photon pile-up percentage of 1.1%, 
1.7%, 3.1%, 4.7%, 6.7%, and 9.2% at 5, 6, 7, 8, 9, and 10 
keV, respectively. Finally, excluding the two brightest columns 
at RAWX= 37 and RAWX= 38, we found a photon pile-up frac- 
tion of 0.2%, 0.6%, 1.6%, 1.9%, 4.8%, and 6.5% at 5, 6, 7, 8, 9, 
and 10 keV, respectively. In this case the 0.5-12 keV Epic-pn 
count rate is 920 c/s. The photon pile-up percentage values with 
the associated errors are reported in Tab. |2] and the plot of the 
fraction patterns of single and double events are shown in Fig. 
|4] After this investigation we choose to extract the Epic-pn spec- 
trum excluding the two brightest columns at RAWX= 37 and 
RAWX= 38 from the source region to avoid any doubt on the 
validity of our analysis. 

Since the high available statistics the EPIC-pn spectrum 
shows evident calibration issues between 2 and 2.5 keV due to 
the Au edge near 2.3 keV, we excluded the EPIC-pn data be- 



4 6 
channel energy (keV) 



10 



Fig. 5. Residuals (in units of cr) in the 0.7-10 keV energy range. 
The adopted continuum is composed by a blackbody plus a mul- 
ticolored disk blackbody both absorbed by neutral matter. In the 
residuals clear emission features are evident at 1, 2.7, 3.35, 3.9 
and in the Fe-Ka region of the spectrum. 



five Gaussian components to fit these emission features. Since 
the line around 2.7 keV is near the hole in the Epic-pn data, we 
fixed the energy of the line at 2.62 keV, imposing that it is as- 
sociated to the Ly-a S xvi transition, furthermore we fixed its 
width at 140 eV. We obtain a x]jd.o.f) = 1.51(269) and a 
improvement of Ax^ - 1775. The significance in units of cr of 
the lines at 1, 2.62, 3.32, 3.9, and 6.7 keV are 14, 7, 15, 9, and 
13, respectively. The best-fit values are reported in Tabs.[3](coI- 
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Table 4. Broad Gaussian Lines between 6 and 7 keV. 



6.5 7 
channel energy (keV) 



Fig. 6. Residuals in the 5.5-7.7 keV energy range. These can be 
fitted adopting a Gaussian profile or, equivalently, a relativistic 
smeared line (diskline in XSPEC). 



umn 2) and |4] (column 2). We identify the broad emission lines 
centered at 1.051!;}™^ keV, keV, and 3.93+^-^^ keV as 

a Fe XXII transition, Ly-ff Ar xvm transition, and resonance tran- 
sition of Caxix, respectively. The centroid of the broader line 
between 6 and 7 keV is not identifiable because we find a line 
centered at 6.80+0.02 keV. We also note that, while the FWHMs 
of the broad lines below 6 keV are between 0.2 and 0.4 keV, the 
broader line in the Fe-Ko' region has a FWHM of 0.6-0.7 keV 
(the line profile is visible in the residuals shown in Fig.|6ll. 

We have also investigated the possibility that the Fe-Ka fea- 
ture is a blending of lines. We fitted the broad Fe-Ka emission 
feature as a blending of Fe xxv (6.7 keV) and Fe xxvi (6.97 keV) 
emission lines adopting two Gaussian profiles and fixing their 
energies at the expected rest-frame values; the values of the fit 
are reported in Tabs.[3](column 3) and|4](column 3). Even in this 
case the Fe xxv and Fe xxvi emission lines appear to be broad, 
having widths of 200 eV (FWHM=0.5 keV). 

Under the hypothesis that the broadening of the lines is pro- 
duced by the relativistic motion of the emitting plasma in the 
accretion disk close to the neutron star, we fit the line profiles 
at 2.62, 3.32, 3.90, and 6.80 keV with relativistic disk line com- 
ponents (DISKLINE in XSPEC). Assuming that the emission 
features are produced in the same region, we impose that the 
inner radius, the outer radius, the emissivity index, and the incli- 
nation angle of the system are the same for all these components. 
Adding these components to the continuum emission we find a 
slightly better fit, ^2^^;(t/.o./.) = L47(268). 

We find that the emissivity index of the accretion disk is 
-2.1, the inner radius is less than 25 gravitational radii (GM/c^, 
hereafter Rg), corresponding to 52 km for a NS mass of 1 .4 Mq, 
and the inclination angle of the system is 41°. The best-fit values 
of the parameters are reported in Tab.|5](column 2); we plot data 
and residuals in the 0.7-10 keV energy band and the coiTespond- 
ing unfolded spectrum in the Fe-Ka region in Fig. [T] and Fig. [8l 
respectively. 

The Fe-Ka line has an energy of 6.76 + 0.02 keV and it is 
not compatible (at 90% confidence level) with the rest-frame 
value of Fexxv transition (6.7 keV). We therefore tried to fit 
the Fe-Ka broad feature adopting two diskline components as- 
sociated to Fexxv and Fexxvi transitions. We fixed the energy 



Parameters 


Single Fe-Ka 


Fe xxv-Fe xxvi 




Em. Line 


Blending 


Efiw (keV) 


6.80 ± 0.02 




CTBroad (sV) 


280!30 




hroad (10-^ ph cm-- s-') 


-0.6 




Equiv. width (eV) 


49+5 




significance (cr) 


13 





^Fexxv (keV) 

o-fe XXV (eV) 

he xxv (10^^ phcm"- s"' 

Equiv. width (eV) 

significance (cr) 

'E'Fexxvi (keV) 

ctfc xxvi (eV) 

he XXVI (10"'' phcm"- s" 

Equiv. width (eV) 

significance (cr) 



6.70 (fixed) 

220!20 
2.6+"-^ 



25 



-0.6 

+ 9 



7 

6.97 (fixed) 
280!"° 
2 0+'-° 

^■"-0.6 

21 ±2 



Note — See Tab. |3] 




4 6 
channel energy (keV) 



10 



Fig. 7. Data, model, and residuals (in units of cr) in the 0.7-10 
keV energy range after adding the disk line components dis- 
cussed in the text. 



values at 6.7 and 6.97 keV, respectively. The best-fit parame- 
ters are reported in Tab. |6] (column 2), in this case we obtain 
xled^d.o.f.) = 1.48(268) and similarparameters of the emissivity 
index and inclination angle of the source. We find an emissivity 
index of -2.2 + 0.2, an inner radius less than 26 Rg, correspond- 
ing to 52 km for a NS mass of 1 .4 M©, and an inclination angle 
of the system of 43+3 

The extrapolated absorbed and unabsorbed flux in the 0.1- 
100 keV energy range was 1.2 x 10"^ and 1.4 x 10"^ erg cm"^ 
s respectively. We obtain an equivalent hydrogen column 
density of 0.74 x 10^^ cm ~^; the same value was obtained by 
IChristian & Swan^ (119971) using the Einstein solid-state spec- 
trometer (SSS; 0.5-4.5 keV), estimating a distance to the source 
of 5 + 1.5 kpc (hereafter we will use this distance). Adopting 
this value as distance to the source the extrapolated unabsorbed 
luminosity in the 0.1-100 keV energy range was 4.3 x 10^^ erg 
s-i. 
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Table 5. Best-fit values of the parameters obtained fitting the 
emission features with relativistic smeared lines (diskline in 
XSPEC). The Fe-Ka feature is modelled with one diskline com- 
ponent. 



Table 6. Best-fit values of the parameters obtained fitting the 
emission features with relativistic smeared lines (diskline in 
XSPEC). The Fe-Ka feature is modelled with two diskhne com- 
ponents. 



Parameters 

N„ (10-2 cm--) 

kTss (keV) 
Nbb(10-') 

kToiSKBB (keV) 

^DISKBB 



Aver, spectrum Spectrum 1+3 Spectrum 2 Parameters 



Aver, spectrum Spectrum 1+3 Spectrum 2 



0.735 ±0.004 0.735 (fixed) 0.735 (fixed) N„ (10^^ cm-^) 



1 797+^016 
i:048!|i 



1 n4Q+''028 

-ft 014 

207^0 



1 780+0 023 

n qoa+oois 

"■='^^-0.017 

i.08!!J:p 

213!« 



1 CI c+0.013 

1.05 ±0.02 

190!^, 



kJsB (keV) 
Nss (10-') 



kT 



N 



'diskbb (keV) 

^DISKBB 



n 7-34+0.002 

'-^ -0.004 



1 Sn 1+0013 

o2» 



1.06 ±0.02 



201 



0.734 (fixed) 0.734 (fixed) 



14 



1 78A+0-025 

-mi 

920+" ''" 
"■=^^"-0.020 

I ns+003 

210+"* 



1 190+0.005 



1 056+"'"' 



186+ ' 



Eifav (keV) 
CTifev/ (eV) 

lu-cV 

Equiv. width (eV) 

P emiss. index 
R,„ (Rg) 

Roi/r (^g) 

9 (degrees) 
Es XVI (keV) 

I5 XVI 

Equiv. width (eV) 
E-Ar XVIII (keV) 



I, 



Mr XVIII 

Equiv. width (eV) 
Ecfl XIX (keV) 



I 



^CaXIX "-"-O? 

Equiv width (eV) 15 ±4 



Epexxv (keV) 

ipeXXV 

Equiv. width (eV) 



1 n'S9 +0 006 

85+^^ 
24+* 
22 ±2 



-2.1 ±0.2 
6.2+"i 



1980'+°^°" 
41.4!^;? 



2.62 (fixed) 
9 9+0.5 
^■^-0.9 
5±3 

3.29 ± 0.03 

11 ±4 

3.89 ± 0.03 

3 6+"^ 
•'•"-0.7 



6.76 ± 0.02 

-'■'-0.6 
61 ±9 

1.47(268) 



1 030+"°"'' 

^•"-'"-0.016 

120+f 
42!f 
32 ±2 



1980 (fixed) 
41.4 (fixed) 

2.62 (fixed) 
3.4+°'' 



,+-J-' 



3.32 ± 0.02 
3.8+°-^ 
12 ±4 



3 89+° "3 

-'•°'-0.04 

3.5 ±0.9 
15 ±5 



6.74 ± 0.04 
5.5 ±0.8 
64± 11 

1.24(271) 



1 062+"°"'* 
66+f 



-2 10+°" 

^•^"-0.16 

10!f 
1980 (fixed) 
41.4 (fixed) 

2.62 (fixed) 
0.9 ± 0.6 

2+3 



-3 90+0.02 
^■^°-Q.p4 
9 7+0.6 
^•'-0.8 

9±4 



3.90 ±0.03 
3.2 ± 0.7 
13 ±5 



6.78 ± 0.03 

c 0+0.7 

59 ± 11 
1.30(271) 



Note — Uncertainties are at the 90% confidence level for a single 
parameter. The parameters are defined as in XSPEC. The continuum 
emission is fitted with a blackbody plus a multicolored disk black- 
body, diskbb. The four broad emission features are modelled adopting 
diskline components. All components are subject to interstellar ab- 
sorption. The intensity of the lines are in units of 10-^ ph cm-^ s-'. 



I 



-ifev (keV) 
"itev (eV) 



17 ± 2 Equiv width (eV) 



y8 emiss. index 

Roi/i (Rg) 
6 (degrees) 

E5 XVI (keV) 

Is XVI 

Equiv. width (eV) 
E-Ar XVIII (keV) 

I^r XVIII 

Equiv. width (eV) 
Ecaxix (keV) 

ICfl XIX 

Equiv. width (eV) 
Ef^xxv (keV) 

If c XXV 

Equiv. width (eV) 

Ef.xxvi (keV) 

If E XXVI 

Equiv. width (eV) 
xUd.o.f.) 



1 r)C9+0.006 
^■"■^ -0.005 

86f 
24+3 
22 ±2 



-2.2 ± 0.2 

9+17 

3600'+'*°°° 



43_3 

2.62 (fixed) 
2.0 ±0.5 
5±3 



a an+0.04 
-'•-"-'-0.03 
T 9+0.5 
-"•^-0.8 

10+"* 



3.896 ± 0.012 

3 4+"-5 

-0.7 

14+^ 



6.70 (fixed) 

4 7+04 

40 ± 10 



6.97 (fixed) 
1 6+"^ 
17+' 

^ '-17 

1.49(268) 



1 r)a9+0.009 

120+f 
41+^ 
31 ±3 



_9 17+0.14 

'-0.19 
6+14 

3600 (fixed) 
43 (fixed) 

2.62 (fixed) 
3.4!1? 
8 ±4 



3.32 ± 0.05 
3.7+?? 



+0.03 
14 



3.89 
3.6+ 

15+ 



6.70 (fixed) 
4.9 ± 1.3 

60!^° 



9+" 
10+2° 



1.25(271) 



Note — See Tab. ID 



1 062+°"°'' 

17.0!?^ 
17 ±2 



-2 10+"'2 
9+20 

3600 (fixed) 
43 (fixed) 

2.62 (fixed) 
0.7 ± 0.6 



3.29 ± 0.04 

2.6!^:? 

3.90 ± 0.03 
3.0 ±0.7 

12+3 



6.70 (fixed) 

3 q+l-O 

40+10 



6.97 (fixed) 6.97 (fixed) 



9 I +0.6 
^■'--o.v 
20 ± 10 



1.30(271) 



2.2. Spectral analysis of the selected time intervals 

To study possible changes of the broad emission line in the Fe- 
Ka region of the XMM spectrum of GX 349+2 we extracted the 
EPIC-pn spectrum from the interval 2 (Spectrum 2) and from 
the intervals 1 and 3 together (Spectrum 1+3). Since the 4-10 
keV lightcurve (Fig. [3j middle panel) indicates a higher count 
rate in the time interval 2 than in the time intervals 1 and 3, 
we expect a harder spectrum coiTesponding to the time interval 
2. The lower statistics in the time selected spectra, due to the 
lower exposure times, requires that we fix some parameters of 
the diskline components. We fixed the outer radius and the in- 
clination angle of the system to the values obtained fitting the 
spectrum corresponding to the whole observation. Initially we 
adopted only a diskline component to fit the Fe-Ko' broad fea- 



ture obtaining an energy value of 6.74+0.04 and 6.78+0.03 keV 
coiTesponding to Spectrum 1+3 and 2, respectively. We note that 
while the line energy is compatible with a Fe xxv transition for 
Spectrum 1+3, it is not for Spectrum 2, that is when the count 
rate in the 4-10 keV energy band increases. The best-fit values 
are reported in Tab. |5] (column 3 and 4). 

As a consequence of this result we fitted the Fe-Ko' broad 
feature in Spectra 2 and 1+3 adopting two diskline components 
with energies fixed at 6.7 and 6.97 keV. Also in this case we 
fixed the outer radius and the inclination angle of the system 
to the values obtained fitting the average spectrum. The best-fit 
values are reported in Tab. |6] (column 3 and 4). We find that the 
main differences between Spectrum 1+3 and Spectrum 2 is the 
value of the blackbody normalization that is larger in Spectrum 
2. Also we note that the flux of the line at 1 keV, the S xvi line, the 
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6.5 7 
channel energy (keV) 



Fig. 8. Unfolded spectrum in the energy range 5.5-7.7 keV; the 
Fe XXV relativistic line profile is shown. 



Ar xvm line, the Ca xix line, and Fe xxv line is lower in Spectrum 
2 than Spectrum 1 +3 while the flux of the Fe xxvi line is larger 
in Spectrum 2 than Spectrum 1+3. In the following section we 
discuss all these results. 



3. Discussion 

We performed a spectral analysis of a 22.5 ks XMM observa- 
tion of GX 349+2 in the 0.7-10 keV energy range. The large 
flux from the source caused a telemetry overflow in the EPIC-pn 
data collection, resulting an eff'ective exposure time of 10 ks. We 
fit the continuum emission adopting a blackbody plus a multi- 
colored disk blackbody and both these components are absorbed 
by neutral matter. Five emission features were clearly visible in 
the spectrum. 

Initially we fitted these features using Gaussian profiles and 
finding their energies at 1.05, 2.62 (fixed), 3.32, 3.9, and 6.8 keV. 
The corresponding widths and equivalent width were 90, 140 
(fixed), 190, 100, and 280 eV, respectively, and 22, 6, 12, 9 and 
49 eV, respectively. We associated the first four lines to L-shell 
Fe xxii- XXIII transition, Ly-a S xvi, Ly-a Ar xvm, and Ca xix res- 
onance transition, respectively. The broader emission feature in 
the Fe-Ko' region was not identifiable, hence we supposed that 
it is a blending of Fe xxv and Fe xxvi emission lines. Fitting the 
broad feature adopting two Gaussian profiles centered at 6.7 and 
6.97 keV we obtained that their corresponding widths were 220 
and 280 eV. We conclude that if the Fe-ka feature is a blending 
of Fe xxv and Fe xxvi emission lines, these lines are intrinsically 
broad. 

BeppoSAX observed three emission lines in th e spectrum 
of GX 349+2 jlaria et al.1 120041 iDi Salvo et aLll2001,) : the first 
one centered at 1.18 + 0.03 keV, interpreted as L-shell Fexxiv 
transition, the second at 2.60+0.06 keV, interpreted as Ly-a S xvi 
transition and, finally, the third at 6.75 ± 0.01 keV, interpreted as 
Fe xxv transition. In this observation we find that the centroid of 
the broad emission line at 1 keV is at a lower energy suggesting 
a L-shell transition of less ionized iron (Fe xxii-xxiii) although 
we cannot exclude a more complex blending of Fe xxi-xxiv. The 
widths and equivalent widths of these three lines are very similar 
comparing our observation and the two BeppoSAX o bservations 
of GX 349+2 (llaria et al.ll2004tlDi Salvo etaPIIOOl in the non- 



flaring state, llaria et akl (|2004|) found, during interval 2 in the 
non-flaring state, widths of < 80, < 120 and 210 ± 80 eV, and 
equivalent widths of 23 + 8, 7^p and 52 ± 1 1 eV corresponding 
to the Fe xxiv, S xvi, and Fexxv emission lines. Furthermore, 
iDi Salvo et al.l (|2001|) obtained equivalent widths for the Fe-Ko- 
line of 71 eV and 34 eV, in the non-flaring and flaring state of 
the source, respectively. From these rough comparisons it seems 
that the source is in a non-flaring state during the observation 
discussed in this work, although we cannot assert this conclusion 
without an accurate analysis of the time variabilities. 

These emission features can be fitted slightly better adopt- 
ing smeared relativistic lines. The large widths of the lines 
suggest that these are produced by reflection of the primary 
spectral component by the accretion disk. Similar broad emis- 
sion features, discussed as smeared Unes were recently ob- 
served in several systems containing a neutr on star, such as 
Ser X-1 (Bhattacharvva & Strohmaver 2007), SAX J1808.4- 
3658 tPapitto et al.. 2009; Cackett et al. 2009a), 4U 170 5-44 
(IDi Salvo et al.ll2005l l2009h. G X 340+0 (iD'Ai et al.ll2009h. 4U 
1636-536 (lPandeletal.1 120081) . X 1624-490 (laria etalj |200^, 
and 4U 1820-30 and GX 349+2 (jCackett et al. 200irand au- 
thors of these papers agree to interpret the broadening as due 
to relativistic and Doppler effects. In the following we show the 
self-consistency of the disk reflection scenario discussing the pa- 
rameters reported in Tab. |6] 

We infer the inner radius of the accretion disk, Rdisk, using 
the normalization value of the multicolored disk blackbody and 
the inclination angle of the source, / = 43^3 deg, obtained by the 
relativistic smeared line, finding that Rdisk is 8 +3 and 8 ±2 km for 
Spectrum 1+3 and 2, respectively. . The blackbody radius, Rbb, 
inferred by the best-fit values of the blackbody normalization is 
4.2 ± 1 .4 and 4.5 + 1 .4 km for Spectrum 1 +3 and 2, respectively, 
too small to be associated to the X-ray emission from the NS 
surface. 

Taking into account possible modifications of the multi- 
colored disk blackbody component due to ele ctron scattering 
(IShakura & Svunvaevlll97l IWhite etaTI Il988l) . flie measured 
color temperature, Tcoi, is related to the effective temperature of 
the inner disk, Tcoi - fTeff, where / is the spectral hardening 
fact or, and Re ff - f^Rmeasured - Adopting / = 1 .7 as estimated 
by Shimura & Takaharal (1 19951) for a luminosity around 10% of 
the Eddington limit, close to the 25% of the Eddington limit esti- 
mated for GX 349+2 in our analysis, the corrected value of Rdisk 
is 24 ± 7 km. The radius of the accretion disk is compatible with 
the value of the inner radius obtained by the relativistic smeared 
lines, that is, assuming a neutron star mass of 1.4 Mq, < 41 km 
and 19^^' km, for Spectrum 1+3 and 2, respectively. In the same 
way, correcting the radius of the blackbody component for the 
electron scattering, we infer that Rbb = 13 ± 4 km. This radius 
is compatible with the NS radius and suggests that the black- 
body component is produced very close the NS surface. In this 
scenario we directly observe the emission from inside the inner 
radius of the accretion disk, that is close the NS surface, and the 
blackbody component is interpreted as the direct emission from 
a compact corona and/or a boundary layer (BL) located between 
the neutron star and the inner radius of the disk. We should ex- 
pect a Comptonized emission from the compact corona and/or 
BL, but if it is optically thick the Comptonized component re- 
sults saturated and mimics a blackbody emission in the adopted 
energy range. 

We observe, for the first time, reflection from the inner accre- 
tion disk that involves not only the K-transitions of iron but also 
ions of lighter elements as Sulfur, Argon and Calcium which all 
correspond to a similar ionization parameter We observe emis- 
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sion lines associated to S xvi, Ar xviii, Ca xix, and Fe xxv, which 
correspond to a ionization parameter Log(f) ^ 3. A similar de- 
tection was recently reported by Di Salvo et al. (2009) analysing 
XMM data of 4U 1705-44. 

Using the relation ^ = L JijieP'), where is the unab- 
sorbed luminosity in the 0.1-100 keV, Ue the electron density of 
the plasma, and r the distance from the source of the emitting 
plasma we can roughly estimate the electron density tie adopting 
as distance from the source the inner and outer accretion disk 
radius, Ri„ and Rout, obtained from the diskline components. We 
obtain that the electron density n,, decreases from 10^2 to 10'^ 
cm"^ going from the inner to the outer radius. These values of 
electron densities ar e in agreement with the values reported by 
IVrtilek et alj ( Il993h . who estimated the electron density of an 



extended cloud to be 10 cm near the equatorial plane at a 
distance from the neutron s tar surface lower than 10^ cm (see 
Fig. 2 in 'VrtileketaIlll993b . 

From the analysis of the two selected spectra we find that 
the extrapolated (not absorbed) flux in the 0.1-100 keV band 
of the multicolored disk blackbody is approximatively constant: 
(0.61 + 0.12) X 10-^ and (0.49 ± 0.05) x lO"** erg cm^^ g-' for 
Spectrum 1+3 and 2, respectively. The blackbody component 
has a larger flux in Spectrum 2 (0.94 ± 0.03 x 10"*^ erg cm"^ s"') 
than in Spectrum 1+3 (0.77+0.02X 10"^ erg cm"^ s"'). This sug- 
gests that the intensity change during the observation is driven by 
the blackbody component. In our scenario the blackbody compo- 
nent is associated to the emission from a BL around the compact 
object and this emission illuminates the surface of the inner ac- 
cretion disk producing the reflection lines observed in the spec- 
tra. We note that the line fluxes are larger in Spectrum 1 +3 than 
in Spectrum 2 except for the flux associated to the Fe xxvi line 
which shows an anticorrelated behaviour. 

Our scenario may easily explain these results: in fact, if the 
emission from the BL increases, the flux illuminating the ac- 
cretion disk increases, then the ionization parameter of the re- 
flecting matter also increases under the hypothesis that the den- 
sity of the matter does not change. This implies that transitions 
of heavier ions are more probable, since lighter ions as S xvi, 
Ar XVIII, Ca xix, and Fe xxv decrease in number while heavier 
ions as Fe xxvi increase in number, giving as results a larger flux 
of the Fe xxvi line and a lower flux of the other lines. 

GX 349+2 is a NS LMXBs belonging to the Z-class sources 
and its behaviour is very similar to the prototype of its class 
Sco X-1. Steeghs & Casares (2002) analysed the optical spec- 
trum of Sco X-1 and obtained a firm mass ratio limit of ^ ^ 0.61 
from the phase-resolved spectroscopy, where q is the ratio of 
the companion star mass to the neutron star mass. If we assume 
that the o ptical per iodicity of 21.9 ± 0.4 hours, observed in GX 
349+2 bv lWachter & Mai-gonid 19961) . is the orbital period of the 
syst em, the mass function of GX 349+2 is /(M) = 0.032 Mq 
(see IWachter & Margon|[T99 ^. We find that the inclination an- 
gle of the system is 43^3 deg. Assuming a neutron star mass 
of 1 .4 M0 we find q ~ 0.56, a value very similar to the upper 
limit estimated for Sco X-1, giving a mass of the companion star 
of M, = 0.78+1}°^ M q typical for a NS LMXB of the Z class 
{q < 0.61 for Sco X-1, [ Steeghs & Ca s"^l2002l: q = 0.34 + 0.04 
for Cyg X-2. ICasares"e t al. 1997). 

Furthermore Sco X-1 shows twin compact radio lobes form- 
ing an angle, with respect th e line of sight, of 44+6 deg (error at 1 
cr, see lFomalont et al.l200ll) . If the jet from the source producing 
the radio lobes is almost perpendicular to the equatorial plane of 
the system then the angle of 44+6 deg corresponds to the inclina- 
tion angle of the system, similar to the inclination angle of 43^3 
deg obtained for GX 349+2 in this work. Our results indicate that 



Sco X-1 and GX 349+2 are ver y similar sources, as already sug - 
gested by other authors (see e.g. lKuuIkers & va"n der Klis|[l998l) . 

Although our scenario looks reasonable, our values of the 
diskline p arameters are not compatible with t he recent ones re- 
ported by ICackett et al.l (|2008|) . ICackett et al.l analysing Suzaku 
data of GX 349+2 fitted the broad emission feature in the Fe- 
Kc region adopting a smeared relativistic line with the outer ra- 
dius fixed at 1000 Rg. They found a line energy associated to 
Fexxvi, an inner radius of 16.5 + 0.8 km, an emissivity index 
of -4.1 + 0.3 and, finally, an inclination angle of the system of 
23° + 1°. The diff'erences in line energy, emissivity index, and 
inner radius of the reflecting skin could be explained suppos- 
ing that GX 349+2 was observed in a different spectral state. 
We find that the Fe-Kor broad emission feature is associated to 
Fe xxv line or to a blending of Fe xxv and Fe xxvi in which the 
Fe xxv line is more intense. Our measurement is similar to the 
previous BeppoSAX o bservations of GX 349+2 (iDi Salvo et al.l 
I2001l:llariaet^l2004 and to the two recent Chandra observa- 
tions of this source dCackett et"aLll2009bl) . In all these observa- 
tions GX 349+2 was in NB/FB and/or in the bottom part of the 
FB. Th e only excep tion is the Suzaku observation of GX 349+2 
(ICackett et al.l 12008 ) where the broad emission line seems as- 
sociated to Fe XXVI suggesti ng that the sou rce was in a different 
state; in fact, as discussed bv lCackett et al.L it is possible that GX 
349+2 was in NB during that observation. 

However, the diff'eren t inclination angles o f the system ob- 
tained in this work and bv lCackett et al.l ( l2008l) cannot be easily 
explained since the inclination angle of the system should not 
depend on the state of the source and should n ot change. Using 
the inclination angle obtained by Cackett et al. we obtain a mass 
ratio q ~ 1.25 and a companion star mass of 1.74^q[^ Mq for 
a neutron star mass of 1.4 M©, that is quite unusual for a NS 
LMXBs. 

Higher statistics spectroscopic studies are needed to put 
tighter constraints on the system parameters of this source, since 
these can give important informations on the whole binary sys- 
tem and on the physics of the accretion flow close to the compact 
object. 



4. Conclusions 

We analysed a XMM observation of GX 349+2 finding the pres- 
ence of prominent emission features associated to S xvi, Ar xviii, 
Ca XIX, Fe xxv and, possibly, Fe xxvi. The emission features can 
be fitted with relativistic smeared lines (diskline). The contin- 
uum emission was fitted adopting a multicolored disk blackbody 
plus a blackbody component. We investigated the scenario in 
which the broad emission lines are formed by reflection from a 
ionized skin in the inner region of the accretion disk illuminated 
by the emission of a compact corona (or a boundary layer) sur- 
rounding the neutron star that we fitted with a blackbody com- 
ponent with a temperature of 1.8 keV. We find a self-consistent 
explanation of our results. The inner radius of the disk black- 
body is 24 + 7 km while the broad lines give an inner disk radius 
less than 40 km from the neutron star center. We discuss that 
the electron density of the reflecting plasma is between 10'^ and 
10^^ cm"^ and that the inclination angle of the source is around 
43°. Finally, we infer that, for an inclination angle of 43°, the 
mass of the companion star is 0.78 M©. 
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